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From a forward-backward path integral, we derive a master equation for the emission and 
absorption of gravitons by a massive quantum object in a heat bath of gravitons. Such an equation 
could describe collapse phenomena of dense stars. We also present a useful approximate Langevin 
equation for such a system. 



I. INTRODUCTION 

Usually, quantum effects are negligible for objects of a large mass which consist of a great number of particles. 
Under special conditions, however, such objects can be of microscopic size. Then the weak gravitational forces may 
sum up coherently. The process of a collapse of a star will be accompanied by an emission of gravitational waves 
which should be treated quantum mechanically in the final stages, e.g. during the formation of primordial black holes. 
In this paper we set up an equation for describing the behavior of a gravitational object of this kind when interacting 
with an equilibrium state of gravitons. 



II. PATH INTEGRAL FOR DENSITY MATRIX 



Consider the Schrodinger equation 



(2.1) 



where H(t) = —h 2 A/2M + U(x) is the Hamilton operator formed from the Laplace-Beltrami operator A = 
g- 1 / 2 (d^gW g 1//2 <9„) containing the inverse g^ v of the metric g^ v {yi,t) and the determinant g =det {g^)- In empty 
space, small fluctuations 2eu fJ , l y(x) of the metric around the Minkowski metric rj^ u describe gravitational waves which 
may be quantized as usual. In the linear approximation, the Einstein action reads 
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where G M1/ is the linearized Einstein tensor H/*" — g^R/2, and e a gravitational coupling constant related to Newton's 
gravitational constant e 2 = 8irG/c 3 . In the radiation gauge with u 0a = 0, the free fields may be expanded into 
creation and annihilation operators for gravitons as 
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where e M;y (k, ±2) are the traceless divergence-free polarization tensors of helicity ±2. The correlation functions at 
temperature T of these fields are: 
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is the projection tensor to the physical polarization states, expressed in terms of transverse projection matrices 
Ajj(k) = Sij—kikj/k 2 of electromagnetism. It is easy to verify the projection property: Atj ) jw(k)A fc ' Tnn (k) = Ay- imn (k). 
As usual, f3 = l/fc^T, and ks is the Boltzmann constant. We have used the four- vector notation x = (x, t) for space 
and time. The operator T is the time- ord ering operator. 

We solve the Schrodinger equation (2.1) with the initial condition (x|-0(O)) = ^(x) and average over the gravitons. 
As a result we obtain a density matrix 



p(x+,x_;i) = (x + |V(t))(V(*)|x_) 
The time evolution of p(x + , x_; t) is described by the evolution equation 



p(x+ a , x- a ; h) = J dx +b dx- b U (x +fc x_ b , i fc |x +a x_ a , t a ) p (x +b , x_ h , t a ) 



(2.6) 
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With the help of the influence functional of Feynman and Vernon (l]-|5| , the kernel of this integral equation is expressed 
by the forward-backward path integral (see Q , and compare with the electromagnetic case in || ) 



U (x +b x^ b ,t b \x +a x- a ,t a ) = I X>x + X>x_exp 

e 2 M 2 r t h 



i 



dt' 



M 



[4-x 2 _] -^(x+) + ^(x_) 



x exp 



2?T 



dt / dt 



e 2 M 2 r t b 



27T 



dt \ dt 



X^_X_^j -^-kl;rnn ^^-f- ^-j- ^ X_ ^ 

i^X^_X_^_ X^_X_^j C~kl:rnn ^*^-f- ^-\- — ^ 



(2.8) 



where x' stands for x(t') and Aki- m n arL d Cki- mn are twice the real and imaginary parts of (2.4). We are assuming the 
wavelengths to be much large than the size of the object, which allows us to neglect the x-dependence in Aki :mn and 
Cki-,mn n- Then we take the angular average 



(Afci ;m „(k)) = Afej ;r 



and find 



and at high temperature: 
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Afei ;mn S(t — t'), 



Introducing the traceless dimcnsionless tensor 
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(2.12) 



the second exponent can be rewritten as 



dt( q ^-q k J)d t (q^+q k _!) R 



dt(q>?-q k J) 2 , 



(2.13) 



where 



e 2 M 2 c 2 k B T 
107m n 



(2.14) 



and the superscript R emphasizes the retardedness of the qki-terms with respect to the g^-terms. This will be 
important to fix the operator order in the master equation to be derived below. Note that 7 is dimensionless, whereas 
w is a frequency. 
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III. MASTER EQUATION FOR DENSITY MATRIX 



Let us derive an operator differential equation for the time evolution of the density matrix. We denote the Hamilton 
operator without gravitational field by 



(3.1) 



and proceeding as in the electromagnetic discussion in Ref. || . The result is the master equation 



' M q kl p~ pq^q 1 " + q™pq"" -q™pq 



kl — 



where x k , x k , q kl , . . . are defined recursively by the commutators x k = i[H, x k ]/Ti, 'x k = i[H,x k ]/h, q 
i[H,q kl ]/h,... . 

This equation will conserve the positivity of the probability as can be assured by bringing it the the extended 
Lindblad form Q 



dtp = - l T[H, P ]-jZY. P + \P%W ~ WW 
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with the Lindblad operators 
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(3.4) 



The master equation (3.2) allows us to calculate the rate of change of an initial eigenstate \i) of the Schrodinger 
equation H\i) — Ei\i) by emission and absorption o f gra vitons to lowest order in G. It is a purely classical rate, 
determined by the matrix elements of the 7-terms in ( |3.2[ ) (see [|| for more details): 

9tPi = -- («|^fci|jf><jf |$ fci l*> — <*|«fci|y> <-f l«fe2]*> — <*lfffc«l*X«l§fcil*> -t-<*l^fcil*> <*l«fc2l*> = -l^Uif{i\m\f){f\qki\i). (3.5) 

The part of the sum with f < i, where the final energy is lower than the initial energy, yields the spontaneous decay 
rate of the system: 



T = 2r/^2wi f {i\q M \f){f\q kl \i) 



(3.6) 



For more details of this identification, in particular the origin of the factor 2, see Ref. g. This decay rate coincides 
with the perturbative quantum-mechanical result obtained from Fermi's golden rule and the interaction operator 
M^ 1 f f * b UkipkPl- 
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In terms of T, the right-hand side of (|3.5|) can be written as 



d t pi = + j^tv t f(i\q k i\f)(f\qki\i) ~ l^uj if (i\q k i\f)(f\qki\i) 

f<i f>i 



(3.8) 



By spontaneous emission, pi decays like e _rt . The rest is due to induced emissions and absorptions contributing to 
the rate of change of pi. 

The total level width is increased further by the induced emissions and absorptions proportional to w in (|3.2|): 
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(3.9) 
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These are due to further induced emission and absorption processes (see Ref. (s) for details) . 

Only the last term in this expressions is caused by the quantum nature of the emitted gravitons. In fact, Eq. (3.7) 
agrees with the classical result for the emitted power by inserting an energy per graviton hajif into the sum and 
replacing ojfj times the square of the matrix elements by the classical expression (M~ 2 dtqki) 2 /2. Recall that according 
to the correspondence principle, we find the quantum-mechanical result for the transition \i) to |/) from the classical 
formula by replacing qij(t) — > (f\qij \i)e~ luJift +c.c. — 2\{f\qij\i)\ cos(iVift — S) from which we deduces the replacement 
of the classical (dtqij) 2 by the temporal average of the quantum-mechanical 2oj 2 j:\(f\qij\i)\ 2 . 

A possible future application of our master equation ( |3.3| ) could be to neutrino oscillation experiments, modifying 
the discussions of Refs. [pU9|. 



IV. LANGEVIN EQUATION 



At high and moderately high temperatures, the fluctuations of y are small due to the last term in the exponent of 
fl2.8p. We therefore introduce average and difference variables of the forward and backward paths in (2.8): 



x = (x++x_)/2, 



y = x+ - x_ 



(4.1) 



and expand V(x±) — V(x± y/2) in y up to the third order. Then the exponent in Eq. (2.8) becomes 
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(4.2) 



where fourth and higher orders in y(t) are neglected. Neglecting also the quadratic terms y(t), extremization with 
respect of y(i) gives the classical equation of motion with radiation damping 



xx 2 H — x (xx) 



= 



(4.3) 



From this equation we can calculate the energy dissipation of the gravitational body due to the coupling to the bath 
of gravitons: 
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Expressed in terms of qu , the right-hand side yields a rate of energy loss 



E 



4G 



Qki- 



(4.5) 



This is the formula derived for classical gravitational radiation. The extra term in the brackets is the gravitational 
analog of Schott's term in the energy balance of electromagnetic radiation damping discussed in [5| JlO| , pT| . 

Interpreting this equation twice over t we can conclude that the time-average of the energy decreases owing to the 
dissipative term on the r.h.s. of (4.4). The quadratic fluctuations of y(£) in the path integral ( |2.10 ) are governed by 
the quadratic exponent (4^2), which can be written as a path integral 



VxDy exp < i 



dt 



M i / 



F(x) 
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where u(t) abbreviates the left-hand side of Eq. (p~3|), 
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The exponent can be completed quadratically to 



(4.7) 
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We now introduce the stochastic variable 



(4.9) 



which has the white-noise correlation function 



{r)i(t)rij{t'))=Mh 2 5 ij 6(t-t') 



5ij6(t - 1'). 



(4.10) 



Ignoring the Jacobian of the transformation x(t) — > u(i), the path integral is seen to be equivalent to the Langevin 
equation 



xx 2 + -x(xx) 2 



dt 



(4.11) 



Note that this Langevin equation is purely classical, since the constants yK and wh? are ^-independent. In electro- 
magnetism, a similar equation has been derived in Ref. |j| solved in Ref. [[12| . 

The omitted Jacobian consists of a product of two Jacobians for the transformations x — > u and u — > r\. For the 
latter, it is equal to Det<9 t f2 1//2 which corresponds to an effective action logDet c^Sl 1 / 2 proportional to (5(0), This 
vanishes in dimensional regularization ji] . The other Jacobian is nontrivial and awaits a proper evaluation along the 
lines of Ref. |l3|], and the result may modify Eq. (4.11). 

The equation is applicable only for high temperatures. If we want to apply it to moderately high temperatures, quan- 
tu m eff ects appear. To lowest order in h, these arise from the first nonleading Taylor expansion term of coth[/3M2k/2] 
in (2.4). This term changes the 5-function in Au-,mn of Eq. (2.11) to 



5{t 



n 2 B 2 - 



(4.12) 



This leads to a change 77 — ^ 77 — h 2 f3 2 if}/2A on the right-hand side of t he Langevin equation ( 4.11 ). 

Let us point out that our derivation of the classical equation (4.4) in the present quantum mechanical context is 
much simpler and more direct than the purely classical derivations in the literature. There one always start from 
extended objects to find the forces exerted upon a body by its own gravitational radiation. The subsequent limit of 
a point particle is quite delicate. For reviews on the subject see Refs. [To|jr^ , p^| 
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